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Abstract 

We present a first-principle study of geometrical and electronic structure of hexagonal single- 
walled silicon nanotubes with a monovacancy or a substitutional defect. The C, Al or P atoms 
are chosen as substitutional impurities. It is found that the defect such as a monovacancy or a 
substitutional impurity results in deformation of the hexagonal single-walled silicon nanotube. In 
both cases, a relatively localized unoccupied state near the Fermi level occurs due to this local 
deformation. The difference in geometrical and electronic properties of different substitutional 
impurities is discussed. 
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Cubic-diamond bulk silicon is a very well-known semiconducting material with an en- 
ergy band gap of 1.2 eV. Unlike the cubic-diamond silicon, it would be hard to form one- 
dimensional single- walled silicon nanotubes (SiNTs) mainly because silicon prefers sp 3 bonds 
to sp 2 bonds. However, Bai and his co-workers l] suggested Si tubes formed by the top- 
to-top stacking of square, pentagonal and hexagonal silicon structures and showed using 
ab initio calculations that the pentagonal and hexagonal SiNTs may be locally stable in 
vacuum and be metallic. Thus, their work provided computational evidences for a possible 
existence of such one-dimensional silicon nanostructures, although those structures have not 
been observed experimentally yet. Therefore, at this point, it is very valuable to theoreti- 
cally investigate the stability and related electronic structure of these silicon nanostructures 
with some defects. 

In this Letter, we report the modification in geometrical and electronic structure of hexag- 
onal single- walled silicon nanotubes with a monovacancy or a substitutional defect. The 
first-principles pseudopotential calculations are carried out based on the density functional 
theory^ within the generalized gradient approximation for the exchange-correlation func- 
tional. The ionic potential is described with the Vanderbilt-type ultrasoft pseudopotential. 
Wave functions are expanded in a plane wave basis set with an energy cutoff of 21 Ry imple- 
mented in the PWSCF code. (4] As a model system, we choose the hexagonal SiNT with a 
monovacancy or a substitutional defect. The C, Al or P atoms are chosen as substitutional 
impurities. All chosen model systems are treated by a supercell with the periodic boundary 
condition. The supercell in the lateral direction is as large as 19.2 A and it has 16 times 
the minimal unit cell (16 x 2.4 A) in the tube axis direction. The structures are relaxed 
until the Hellmann-Feynman forces are smaller than 0.05 eV/A. The Brillouin-zone integra- 
tion is done within the Monkhorst-Pack scheme using 1x1x4 k-point sampling. To 
understand the features of charge distribution of our models, we perform the Mulliken pop- 
ulation analysis using a numerical atomic orbital basis set in the OpenMX codefl Q with 
a kinetic energy cutoff of 150 Ry. For this Mulliken population analysis, norm- conserving 
Kleinman-Bylander pseudopotentialsfsl, [9J are employed. 

First, we calculate hexagonal SiNTs with or without a monovacancy shown in Fig. [U If a 
single Si atom is removed from the nanotube, two of four Si atoms around the vacancy rebond 
(bond length of 2.42 A) to form a pentagonal cross section. In the case of cubic-diamond 
bulk silicon, the lattice constant is 5.43 A, which corresponds to the Si-Si bond length of 2.35 



A. Note that a buckled honeycomb planar structure [k| of silicon has the average Si-Si bond 
length of 2.2 A. The rebonded two atoms are perpendicular to the tube axis. Although other 
two Si atoms have dangling bonds even after the relaxation, their dangling bond characters 
become weak since the two Si atoms make strong bonds with three adjacent Si atoms (bond 
length of 2.32-2.36 A). The formation energy (= E perfect - E vacant - E Si atom) of a single 
vacancy in an SiNT is 16.3 eV (4.07 eV per bond). In the optimized structure, the Si atom 
on the opposite side to the monovacancy is protruded by ~1 A from the surface of the 
nanotube. The bonding angle with respect to this protruded Si atom at the pentagonal 
defect is 89.6° as shown in Fig. 1(c). Molecular dynamics (MD) simulations revealed that 
pentagonal and hexagonal SiNTs can remain stable. [lo| Our MD simulation result shows 
that the SiNT is somewhat distorted but not broken when a monovacancy is formed at 600 
K. 

Such locally large distortion may affect the band structure and charge distribution of the 
SiNT. In terms of the Mulliken population analysis, we find an interesting feature. The 
protruded Si atom on the opposite side to the vacancy obtains 0.12 e and its adjacent Si 
atoms lose electrons: Each of two adjacent Si atoms in the tube axis direction loses 0.025 e, 
while each of two in the circumference direction loses 0.035 e. On the other hand, the two 
rebonded Si atoms lose 0.073 e, respectively. 

As mentioned above, Fig. [2(a) clearly shows that the hexagonal SiNT has the metallic 
character. |j| When a monovacancy is formed in the SiNT, the energy band structure is 
changed accordingly. There are prominent localized states above the Fermi level. Especially, 
a localized state (state A) originating from the structure distortion gives a flat band around 
0.2 eV above the Fermi level as denoted by a downward arrow in Fig. [2(b). The isodensity 
surface plot of this unoccupied localized state (ULS) is presented in Fig. [3](a). In addition, a 
pentagon as a topological defect also gives rise to two ULSs (states B and C) near +0.6 eV 
as depicted by upward arrows in Fig. 0(b). Interestingly, Figs. [3(b) and (c) show that the 
localized states by a pentagonal defect are split into even and odd parity states with respect 
to the mirror planes (a c and a n ) containing and normal to the nanotube axis, respectively. 

Next, we present the geometrical and electronic structure of a hexagonal SiNT with a 
carbon, aluminum or phosphorus as a substitutional impurity, shown in Fig. [H In the case 
of the C doping, the C-Si bond length (Cl-Si2 and Cl-Si3) parallel to the tube axis is 1.93 
A and that (Cl-Si4 and Cl-Si5) perpendicular to the axis is 1.92 A. The C atom is buckled 



inward by ~0.4 A and the Si8 atom on the opposite side to the C dopant is protruded 
outward by ~0.6 A. Owing to the strain effect, Si-Si bond breaking occurs near the defect 
site. The bond angle between C and adjacent Si atoms, ZSi2ClSi3, parallel to the tube axis 
are 165.8° and ZSi4ClSi5 perpendicular to the axis is 120.6°. For the Al doping, the Al-Si 
bond lengths parallel and perpendicular to the tube axis are 2.58 and 2.61 A, respectively. 
Here the Al atom is buckled inward by ~0.7 A and the Si atom on the opposite side to 
the Al dopant is buckled inward by ~0.1 A. The bond angle between Al and adjacent Si 
atoms parallel to the tube axis are 133.1° and that perpendicular to the axis is 140.3°. In 
contrast, the P dopant shows an outward buckling by ~0.3 A. Such a trend is shown for 



boron nitride nanotubes. 



Ill ] where nitrogen atoms are buckled outward and boron atoms 



are inward, which is associated with unpaired electrons. For the P doping, the Si atom on 
the opposite side to the P dopant is buckled outward by ~0.3 A. The P-Si bond lengths 
parallel and perpendicular to the tube axis are 2.50 and 2.30 A, respectively. The bond angle 
between P and adjacent Si atoms parallel to the tube axis are 145.5° and that perpendicular 
to the axis is 102.3°. Each impurity (C, Al or P) atom results in more or less flat bands, 
i.e., localized states, (not shown here) about 1.5 eV below or above the Fermi level. Note 
that a boron or nitrogen atom doped in the carbon nanotube also gives the same feature. 
When substituted, the C atom donates 0.25 e to the SiNT. The protruded Si atom on the 
opposite side to the C atom obtains 0.09 e. In contrast, the P atom obtains 0.23 e from the 
SiNT and the protruded Si atom on the opposite side to the P atom obtains 0.014 e. In the 
case of the Al impurity, there is no practical electron transfer between the SiNT and the Al 
atom (0.003 e). It is thus found that the large displacement of the Si atom opposite to a 
substitutional impurity is associated with the large electron transfer. 

Interestingly, Fig. shows the relatively flat bands around 0.1 — 0.2 eV above the Fermi 
level for the substitution of C, Al and P, respectively, which are denoted by red arrows. We 
find that these localized states originate from the structure distortion, so the origin of these 
states is the same as that of the flat band (denoted by a downward red arrow in Fig. [2](c)) 
of ULS obtained by the presence of vacancy. Especially, since the geometrical structure of 
SiNT with C-doping in Fig. |4] is almost the same as that of SiNT with vacancy in Fig. [T](c) 
except the substituted carbon atom itself, the isodensity surface plot of the ULS in C-doped 
SiNT is expected to be almost the same as that shown in Fig. [3J 

In summary, structural and electronic property changes caused by various vacancy or 



substitutional defects are studied in the hexagonal SiNTs using ab initio pseudopotential 
calculations. The C, Al or P atoms are chosen as substitutional impurities. It is found that 
the defect such as a monovacancy or a substitutional impurity results in deformation of the 
nanotube. The relatively localized states near the Fermi level shown in both cases with a 
monovacancy or an impurity occur due to this local deformation. When we consider the 
realistic situation of the randomly distributed defects, there would be somewhat broadened 
density of states originating from almost degenerate localized states. 
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FIG. 1: Schematic ball-and-stick models of (a) a pristine hexagonal SiNT, (b) an SiNT with a 
monovacancy before relaxation, and (c) an SiNT with a monovacancy after relaxation, respectively. 




FIG. 2: Band structures for the SiNT structures corresponding to (a) a pristine SiNT and (b) an 
SiNT with a monovacany after relaxation. States A, B, and C are ULSs. The Fermi level is set to 
zero. Localized states occur originating from a monovacancy. 




state A states state C 



FIG. 3: Isodensity surface plots of the ULSs at T that comes from the distortion due to the presence 
of vacancy: Top and side views for a ULS (state A) near +0.2 eV due to the local deformation are 
shown in (a), and those for two ULSs (states B and C) near +0.6 eV due to the pentagonal defect 
are in (b) and (c). The values for the red and blue isodensity surfaces are ±0.02 e/ajj, where the 
sign is that of the wave function and ao = 0.529 A, the Bohr radius. 
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FIG. 4: Schematic ball-and-stick models with a substitutional impurity in a hexagonal SiNT. The 
grey, black, pink, and yellow balls represent a silicon, carbon, Al, and P atoms, respectively. 
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FIG. 5: Band structures for carbon, aluminum and phosphorus atom substitutions. Arrows indicate 
ULS originating from the deformation near the Fermi level. 



